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Although many of works have been done in multivalued quantum logic synthesis, the question whether mul¬ 
tivalued quantum circuits are more efficient than the conventional binary quantum circuits is still open. In 
this article we devote to the optimization of generic multivalued quantum circuits. The multivalued quantum 
Shannon decompositions (QSD) are improved so that the circuits obtained are asymptotically optimal for all 
dimensionality d. The syntheses of uniformly multifold controlled R y rotations are also optimized to make the 
circuits further simplified. Moreover, the theoretical lower bound of complexity for multivalued quantum cir¬ 
cuits is investigated, and a quantity known as efficiency index is proposed to evaluate the efficiency of synthesis 
of various quantum circuits. The algorithm for qudit circuits given here is an efficient synthesis routine which 
produces best known results for all dimensionality d, and for both cases the number of qudit n is small and that 
is asymptotic. The multivalued quantum circuits are indeed more efficient than the binary quantum circuits. 
The facts, the leading factor of the lower bound of complexity for qudit circuits is small by a factor of d - 1 
in comparison to that for qubit circuits and the asymptotic efficiency index is increased with the increase of 
dimensionality d, reveal the potential advantage of qudit circuits over generic qubit circuits. The generic n- 
qudit circuits with d > 5 and generic two-ququart circuits synthesized by the algorithm given here are practical 
circuits which are more efficient than the most efficient qubit circuits. 

PACS numbers: 03.67.Lx. 03.67.Ac 


I. INTRODUCTION 

Enormous progress has been made in the held of quan¬ 
tum information science over the past two and a half decades. 
Most approaches to quantum information processing use two- 
level quantum systems (qubits). However there is increasing 
interest in exploiting protocol with multilevel quantum sys¬ 
tems (qudits) [ l ~9i]. The simplest multilevel system, the three- 
level quantum system, is called a qutrit, the four-level quan¬ 
tum system is called a ququart. The multivalued quantum in¬ 
formation is exciting because quantum systems usually have 
multi-levels, it enables us to full use of various resources. 

In quantum computing, the algorithms are commonly de¬ 
scribed by the quantum circuit model. The process of con¬ 
structing quantum circuits by some elementary components 
is called synthesis. The complexity of quantum circuit can 
be measured in terms of the number of elementary gates re¬ 
quired. A large amount of work in these areas has been done 
for binary quantum computing fl^NETIl . The CNOT gate is 
one of most widely used two-qubit elementary gate. It has 
been shown that the CNOT gate with one-qubit gate is univer¬ 
sal for qubit quantum circuits [illllll]. The best result so far 
for the synthesis of generic qubit quantum circuits are given 
by Shende et al. based on quantum Shannon decomposition 
(QSD) |d. 

Although many of works also have been done in multi¬ 
valued quantum logic synthesis ll22l - l29ll . the works on this 
area are still far from complete. The results obtained cannot 
show the advantage of using multilevel quantum systems in 
the complexity of quantum logic synthesis. Which gate is cho¬ 
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sen as the two-qudit elementary gate of the qudit quantum cir¬ 
cuit is a crucial issue for multivalued quantum computing, and 
there have been many proposals. In our recent previous work, 
the generalized controlled X (GCX) gate has been proposed 
as the two-qudit elementary gate for the multivalued circuits 
MM- We generalize QSD, the most powerful technique for 
the synthesis of generic qubit circuits, to the multivalued case. 
Based on the GCX gate, using the multivalued QSD, we ob¬ 
viously improve the results of the synthesis of qudit quantum 
circuit & But there are still some problems. One is that the 
quantum circuits built by multivalued QSD algorithm are not 
asymptotically optimal except for the dimensionality of qudit 
d is a power of two. It was not clear whether we can build 
efficient quantum circuits for the qudit d is not a power of two 
as that d is a power of two. The other is that the multivalued 
quantum circuits in Ref. 0 do not show obvious advantage 
over the circuits for binary systems. The problem whether the 
multivalued quantum circuits can be more efficient than the 
binary circuit is still open. 

In this article, we devote to optimizing multivalued quan¬ 
tum circuits and to solving the problems stated above. The 
multivalued QSD for the qudit d is not a power of two is opti¬ 
mized so that the synthesis of quantum circuits for these qudit 
also is asymptotically optimal. The synthesis of the uniformly 
multifold controlled R y rotations is also optimized to make 
the circuits further simplified. The theoretical lower bound 
of complexity for qudit quantum circuits is investigated. A 
quantity known as efficiency index is proposed to evaluate the 
efficiency of synthesis of generic n-qudit circuits. The results 
and comparison show that algorithm given here is most effi¬ 
cient qudit synthesis routine so far which produces best known 
results in all respects. The multivalued quantum circuits are 
indeed more efficient than the binary quantum circuits. 

The article is organized as follows: The lower bound of 
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complexity for qudit circuits is investigated in Sec. M The 
leading factor of the lower bound of complexity for qudit cir¬ 
cuits is small by a factor of d— 1 in comparison to that for qubit 
circuits. The optimization of the multivalued QSD and uni¬ 
formly multifold controlled R y rotations, the structure and the 
GCX gate count of optimal qudit circuits are given in Sec. UHl 
The efficiency of synthesis of quantum circuits is discussed in 
Sec. UV] The quantity in term of efficiency index is proposed 
in this section. The asymptotic efficiency index is increased 
with the increase of dimensionality d for these circuits. The 
efficiency indexes of generic n-qudit circuits with d > 5 and 
generic two-ququart circuits given here are higher than that 
of the most efficient qubit circuits. Finally, a brief conclusion 
and future work are given in Sec. [V] 


II. LOWER BOUNDS OF GCX GATES 


The theoretical lower bound of complexity for generic qudit 
circuits is (d 2n -n(d 2 - 1)- l)/(4 (d- 1)). The lower bound of 
complexity for generic qubit circuits is (4" - 3 n - l)/4 0. 
The leading factor of the lower bound of complexity for qudit 
circuits is small by a factor of d - 1 in comparison to that for 
qubit circuits reveals the potential advantage of qudit circuits 
over the qubit circuits. 

III. OPTIMAL SYNTHESIS OF MULTIVALUED 
QUANTUM CIRCUITS 

A. Synthesis Based on Multivalued QSD jU] 

The hi st phase of multivalued QSD is to use the cosine-sine 
decomposition (CSD) |32j]. Let the m x m unitary matrix W 
be partitioned in 2 x 2 block form as 


The GCX gate (denoted as GC X(m —> X (lJ> )) is a 

controlled- U two-qudit gate which implements the X ( ' J> op¬ 
eration on the target qudit iff the control qudit is in the state 
I m), where X™ = |*>0‘l + \j)d\ + 2 ] k ^j \k)(k\. The GCX gate 
essentially is a CNOT gate. For a multilevel quantum system 
which forms a qudit, two levels in the system forms a qubit. 
If a two-qubit CNOT gate is realized in two such systems, a 
GCX gate is naturally obtained. The number of the GCX gates 
required can be used as a unified measure for the complexity 
of various quantum circuits llTlIl . 

In quantum computing, the quantum circuit is a unitary 
transformation on the quantum states. A generic n -qudit quan¬ 
tum circuit is fully determined by d 2n - 1 real parameters (up to 
a phase factor). Here the qudit circuits are constructed by us¬ 
ing GCX gates and arbitrary one-qudit gates. The GCX gates 
do not introduce any parameters, but they are a kind of bar¬ 
riers that separate one-qudit gates so that they cannot merge 
into a resulting one-qudit gate for each qudit. Intuitively, ev¬ 
ery GCX gate can be accompanied two one-qudit gates, one 
is for control qudit G\ and the other for target qudit G 2 , ap¬ 
plied after every GCX gate. Two one-qudit gates can carry 
2 (d 2 - 1) real parameters, but only part of them can be carried 
by a GCX gate. 

The one-qudit gate corresponds to a SU(d) group, and the 
d 2 - 1 parameters correspond to the bases of suid) algebra. 
Without loss of generality, we consider the GC X{d - 1 —> 
X(d-2,d-i)'j g ate anc j use qjg na m ra l bases of u{d) algebra \i)(j\, 
where i, j e 0,1, • • • ,d- 1. There are 2 (d - 1) bases which do 
not commute with the GCX gate in G 1 , they are | d - 1 )(i\ and 
| i)(d — 11, where i e 0,1, ■ • • ,d — 2, so the gate can separate 
2 (d - 1) parameters in G\. There are Aid - 1) such bases in 
Gi, they are | d - l)(i|, | d - 2){i\, where i e 0,1, • ■ ■ ,d — 1 and 
| j)(d — 11, | j)(d - 2|, where i € 0,1, • • ■, d - 3. But there are 
2 (d - 1) linear combinations of them which commute with the 
GCX gate, they are \d-2)(i\+\d-\)(i\, where i € 0,1, • • • ,d-l, 
and | j)(d -21 + \j)(d - 11, where j € 0,1, • • •, d - 3. The GCX 
gate also can separate 2 (d - 1) independent parameters in G 2 . 
Each GCX gate can bring at most 4 (d - 1) parameters. For a 
generic «-qudit quantum circuits, there are d 2 " - n(d 2 — 1) — 
1 parameters which need to be brought by the GCX gates. 
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with 2 r < m. Here r is called as partition size. Let W = UYV 
be CSD of the matrix, then 



r m - r 



( 2 ) 


r r 
r (C -S 
r =r sc 

m — 2 r ,0 0 


in — 2 r 
0 ' 
0 


(3) 


r m - r 


V = 


r 

m - 


r 




(4) 


where C and S are diagonal matrices of the 
forms C = diag {cos 6 \, cos 62, ■ • •, cos 6 ,•} and 
S - diag {sin 8 \, sin 82, ■ ■ ■, sin (9, ), I is the (m - 2r) x (in - 2r) 
identity matrix, and I is called cosine-sine matrix. An n-qudit 
gate corresponds to a d n x d" unitary matrix. The synthesis 
of qudit quantum circuits based on CSD was first proposed 
by Khan et al. | {27ll28ll . There, they choose the partition size 
r = d "~ 1 at each recursion level. Different from the Khan et 
fl/.’s method, we choose the partition size r = [dl2\d n ~ x for 
the first level decomposition, then r = [d/4]d" 1 for the sec¬ 
ond level decomposition, and r = [d/2 k ]d" for the kd th level 
decomposition, here [a\ denotes the integer part of a. After k 
levels (log 7 d < k < log 2 d + 1) of decomposition, d n ~ l x d n ~ l 
block diagonal matrices are obtained. The block diagonal 
matrices correspond to uniformly controlled (n - l)-qudit 
(uf\\(U n 1 )) gates; the cosine-sine matrices corresponds to 
uniformly (n - l)-fold controlled R y (uA„-i(R v )) rotations. 

The second phase of multivalued QSD is the further decom¬ 
position for the uniformly controlled (n - l)-qudit gate. It can 




3 


be decomposed into d copies of in - l)-qudit gates and d - 1 
copies of controlled (n - l)-qudit diagonal (Ai(A' !_1 )) gates. 
In qubit case, the uniformly controlled in - l)-qudit gate is 
decomposed into a pair of (n - l)-qudit gate and a Ai(A' !_1 ) 
gate, and it is equivalent to the decomposition of the block di¬ 
agonal matrix in QSD. So the decomposition given here is a 
generalization of QSD for qubit case. 

The synthesis of a generic n-qudit gate involves three kinds 
of component: (n - l)-qudit gates, Ai(A" _1 ) gates, uA n -\(R y ) 
rotations. The (n - l)-qudit gates can be further decomposed 
in similar ways. So we can construct a generic n-qudit quan¬ 
tum circuit by a recursive way. All the component elements 
required can be efficiently synthesized based on GCX gate. 


B. Optimization of the Synthesis Stated Above 


Optimizing the decomposition of matrices: When d is not 
a power of 2, there are identity submatrices in the cosine-sine 
matrices of CSD. We can rearrange and change the block di¬ 
agonal matrices of CSD to reduce the numbers of two com¬ 
ponents, the (n - l)-qudit gate and the A|(A" ') gate. The 
number of (n - 1 )-qudit gates can be reduced to the minimum 
d 2 . 

For example, an n-qutrit gate, corresponding to a 3" x 3" 
unitary matrix, can be decomposed as follows 

W X =AY X BT Q CT 2 D. (5) 

with 
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where each block matrix in the decomposition above is of size 
3' 1 * 1 x 3" 1 . We can rewrite it as that 


Wi = A'Y\B'Y 0 C'Y 2 D'. (8) 

with 
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where t/' = U~ l U 3 , X\ = X~ l U~ l Ui, V' = V~ l X- x X 3 V 3 , 
Y\ = Vd l X\. For matrices A, B, C, D and D', each of them 
corresponds to two AAA"' 1 ) gates and three in - l)-qutrit 
gates; for A', B', C' each matrix corresponds to one Ai(A" _1 ) 
gate and two in - l)-qutrit gates. The optimal synthesis of a 
generic n-qutrit circuit gate involves nine (n - l)-qutrit gates 
and five A|(A" *) gates, three in - l)-qutrit gates and three 
A (A" *) gates less than that in the original synthesis. 

For example again, taking d - 6, the n-qudit circuit, corre¬ 
sponding to a 6" x 6" unitary matrix, can be decomposed as 
follows 


W 2 = AY 1 BY 2 CY 3 DY 0 EY 4 FY 5 GY 6 H (10) 
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where each bock matrix in Eq. CD is of size c f- x and 
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A = diag\A ,, A 2 , A 3 , A 4 , A 5 , A 6 ), 

B = diag{I , Z?2 9 B 3> £ £5, B 5), 

C = diag{C\ C2, C3, C4, C5, C6), 

D = diag{I , D 2 , £>3, £ £>5, £>6}, ( 13 ) 

where each block matrix in Eas. (112m~3l i is of size 6" _1 x 6' !_1 . 
The second half of expression in Eq. (ITOl) . TT4FF5G16//. has 
same form as the first half, ,41) BY 2 CY 2 D. The first half can 
be rewritten as 

AriBr 2 cr 3 D = aTiBT 2 ct 3 d' (14) 

with 

A' = ( I®A 2 )diag{I,I,A' 3 ,I,A' 5 ,A ' 6 ), 

B' = (I® B 2 )diag{B\,I , /, B\, B' 5 , /), 

C’ - (I® C 2 )diag{l, /, C 3 , /, C' 5 , Cg), 

D' = diag{D\,D 2 ,D 3 ,D\,D 3 ,D(,}. 


( 15 ) 
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Here A' = A"'A 3 , A' = Aj'Ag, A^ = A^A 6 , B\ = B^A~ l A u 

b' 4 = b^a- 1 a a , b' 5 = b- 1 b 5 , q = c~ 1 b~ 1 b 3 c 3 , q = 
q'c 5 , q = q^j^eCe, q = q'Ci and q = q‘c 4 . 

The second half of the expression can be processed in same 
way. The optimal synthesis of a generic n-qudit circuit gate 
with d - 6 involves 36 (n - l)-qudit gates and 28 A^A"" 1 ) 
gates, eight (n - l)-qudit gates and eight Ai(A" _1 ) gates less 
than that in the original synthesis. 


TABLE I: Numbers of three components for optimal n-qudit quan¬ 
tum circuits. 

d 

component 3 4 5 6 7 8 

(n - l)-qudit gate 9 16 25 36 49 64 

AAA"- 1 ) 5 12 17 28 41 56 

uK-dRy )3 6 10 15 21 28 


The numbers of three components needed to construct op¬ 
timal generic multivalued quantum circuits are listed in Tab. 
U A cosine-sine matrix for a qudit system can involve several 
sets of uA k (R y ) rotation. To reduce the number of (n — 1)- 
qudit gates to its minimum d 2 is essential for the asymptotic 
optimality of the synthesis. The number of the GCX gate in 
these d 2 ( n - l)-qudit gate components account for vast ma¬ 
jority of the GCX gate count of an n-qudit gate if n is large. 
For example£the number of the GCX gate in nine 2-qutrit gate 
components account for 68.02% of the GCX gate count of a 
generic 3-qutrit circuit. Whereas the number of GCX gate in 
nine 7-qutrit gate components account for 99.87% of the count 
of a generic 8-qutrit circuit. The numbers of two other compo¬ 
nents may be neglected if n is enough large. Hence the synthe¬ 
sis obtained here is asymptotically optimal, which means that 
the generic n-qudit circuit can be synthesized asymptotically 
by 0(ad 2 ") two-qudit elementary gates, here a is a constant. 

Optimizing the uniformly multifold controlled R y rotations: 
The optimization of the qubit uA n -\(R y ) rotations has been 
given in Ref. [19] by using controlled-Z (CZ) gates. To opti¬ 
mize the qudit u/\ n -\(R y ) rotations needs the high dimensional 
counterpart of CZ gate. The multivalued extension of Z oper¬ 
ation is a one-qudit operation Z [mI which is specified by that 
Zl™] = Y,k*m - \m){m\. There are only d different forms 
of Z 1 "' 1 operation for a qudit (d = 0,1, • • •, d - 1, respectively), 
whereas there are d(d - l)/2 forms of X <IJ> operation for the 
qudit. Like the GCX gate, the generalized controlled Z (GCZ) 
gate (denoted as GCZ (m - m')) is defined as a controlled 2- 
qudit gate which implements the Z 1 '" 1 operation on the target 
qudit iff the control qudit is in the state | m). It is specified by 
GCZ(m - m') = 2, 7?w | ij)(ij\ - \mm')(mm'\. Like the qubit 
case, the control qudit and target qudit of the GCZ gate are 
changeable, and using two generalized Hadamard gates, the 
GCX and GCZ gates can be transformed each other. The gen¬ 
eralized Hadamard gate is a one-qudit gate specified by 
that H<® - £ k +ij\k)(k\ + (| i)(i\ + |i)01 + \j)(i\ - |7)0I)/V2. 
The circuit representation of GCZ gate and its transformation 
relation with GCX gate are shown in Fig. [T] 

The statements and Fig. 14 in Appendix C of Ref. iH 
still hold for qudit «A„_i(R v ) rotations if all GCX(m —> X <lJ> ) 



FIG. 1: GCZ gate (a) and transformation between GCX gate and 
GCZ gate (b). 


gates are replaced with GCZ (m - j) gates. Thus a set of qu¬ 
dit t<A„_i(R v ) rotation may be implemented with 2d" 2 (d - 1) 
GCZ gates, of which d -1 GCZ (m — j) gates (m = 1,2, • • ■, d- 
1, respectively) may be moved furthest to the right (or the 
left). The rightmost d - 1 GCZ gates produce a diagonal 
gate which may be absorb into the neighboring uniformly con¬ 
trolled (n -1 )-qudit gate. The reason why the d- 1 GCZ gates 
are able to be moved furthest to the right whereas only one CZ 
gate is able to move like that in qubit case is that the controlled 
gates with different control basis states can exchange one an¬ 
other. This saves d - 1 two-qudit elementary gates for each set 
of qudit uX k {Ry) rotation, totally save (d 2{n ~ V) - 1 )n u x/(d + 1) 
GCX gates for a generic n-qudit circuit, where n„R is the num¬ 
ber of uA n -\(Ry) component in the circuit and is given in Tab. 
Q] In the practical process of optimization, it should optimize 
f/A„ i(Rf) rotations first, then to optimize the decomposition 
of matrices. 


C. Structure and GCX Gate Count of Optimal Circuits 

The optimal quantum circuit of generic n-qudit circuits in¬ 
volves d 2 (n - l)-qudit gates, which are separated by A^A' 1 " 1 ) 
gates or circuits for cosine-sine matrix. It involves d 2 - 2“ 
AAA"" 1 ) gates and 2 K - 1 circuits for cosine-sine matrix. In 
multivalued case, a circuit for cosine-sine matrix usually in¬ 
volves several sets of uK n -\{R y ) rotations. The structure of a 
generic n-qutrit circuit is illustrated in Fig. [2] 

The GCX gate count of the optimal multivalued quantum 
circuits given here is tabulated in Tab. |TT] F°r comparison the 
count before the optimization is also given. The results are ob¬ 
viously improved after optimization, especially for the case d 
is not a power of two. The circuits optimized have asymptotic 
optimal features for all dimensionality d, whereas the circuits 
before optimization are not asymptotically optimal except for 
d is a power of two. 


IV. EFFICIENCY OF SYNTHESIS OF QUANTUM 
CIRCUITS 


To evaluate the efficiency of synthesis of generic n-qudit 
circuits based on GCX gates, we propose a quantity known as 
efficiency index ( Sf ) which is defined by Sf — d 2n /N n , where 
N„ is the number of GCX gate required to synthesize the n- 
qudit circuit. The quantityis the average number of param- 
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FIG. 2: Structure of a generic ;i-qutrit circuit. Here the small square (□)denotes the uniform control, the slash (/) represents multiple qutrits 
on the line. 


TABLE II: The GCX gate count for the synthesis of qudit quantum circuits obtained using the multivalued QSD. In each cell, the upper line 
denotes the count before optimization 13 ill , the bottom line denotes the count optimized. 


d 


3 

4 

5 

6 

7 

8 

44 

108 

272 

510 

828 

1176 

26 

90 

176 

355 

618 

980 

692 

2232 

10256 

25860 

52740 

85456 

344 

1926 

5216 

15565 

53856 

72716 

6860 

37800 

336144 

1158720 

2965788 

5551504 

3458 

32886 

136576 

577705 

1797210 

4735948 

83924 

613248 

10796560 

51109320 

166400964 

355955600 

32028 

534582 

3445576 

20902225 

88346400 

303759820 

1011932 

9854970 

345689872 

2.25 x 10 9 

9.32 x 10 9 

2.28 x 10 10 

291638 

8587062 

86295576 

752674425 

4.33 x 10 9 

1.95 x 10 10 

12157748 

1.58 x 10 s 

1.11 x 10 10 

9.90 x 10 10 

5.22 x 10“ 

1.46 x 10 12 

2634932 

1.35 x 10 s 

2.16 x 10 9 

2.71 x 10 10 

2.12 x 10“ 

1.25 x 10 12 

1.46 x 10 s 

2.52 x 10 9 

3.55 x 10“ 

4.36 x 10 12 

2.92 x 10 13 

9.34 x 10 13 

23744984 

2.20 x 10 9 

5.39 x 10 10 

9.75 x 10“ 

1.04 x 10 13 

8.00 x 10 13 


eters carried by each GCX gate. The larger the quantity M, the 
more efficient the synthesis of the circuit is. For the quantum 
circuit synthesis which is asymptotically optimal, there is an 
asymptotic efficiency index .:if asy which is the efficiency index 
when n is enough large. The asymptotic efficiency indexes for 
the optimal synthesis given here are listed in Table III. From 
the Table, it can be seen that the .^ asy is increased with the 
increase of dimensionality d. 

TABLE III: The asymptotic efficiency indexes for the optimal syn¬ 
thesis of multivalued quantum circuits. 


d 

3 

4 

5 

6 

7 

8 

^asy 

1.81 

1.95 

2.48 

2.89 

3.19 

3.53 


There are several previous works on the synthesis of multi¬ 
valued quantum circuits based on elementary gates. Based on 
the controlled increment (CINC) gate, the synthesis by using 
the spectrum decomposition algorithm is investigated in Refs. 
1123112511 . It is asymptotic optimal, which has leading factor 2 
for the CINC account of the synthesis. Using the GCX gate 
as the two-qudt elementary gate instead of the CINC gate, the 
synthesis is greatly simplified ll3ltl . The synthesis simplified 
still has leading factor 2 but for the GCX account. So it’s 
/4? asy is equal to 0.5 for all dimensionality d, which less than 
all values of /^ asy in Tab. [Ill] 

Based on CINC gate, the syntheses by using the CSD with 


balanced partition are investigated in Ref. [28!]. The circuits 
synthesized by this method are simpler than those by using the 
spectrum decomposition if n is small, but they are not asymp¬ 
totically optimal except for d is a power of two. The results 
of this work are given in Tab. HV] It needs d - 1 GCX gates 
to synthesize a CINC gate EOT Comparing the data in Tab. 
[II]and that in Tab. Uv] considering the CINC gate itself has 
complex construct, it can be seen that the synthesis of quan¬ 
tum circuits given in this article are much more efficient than 
that in Ref. [28] even if the n is very small. 

Li et al. propose a two-ququart gate, termed the controlled 
double-not (CDNOT) gate, for four level quantum systems. 
Based on the CDNOT gate, they investigate the synthesis of 
ququart quantum circuits by using QSD method [29], the re¬ 
sults are tabulated in Tab. [V] A CDNOT gate is a two-ququart 
controlled gate which implements the cr v ® L operation on 
the target ququart iff the control ququart is in the state | m), 
m e 2, 3, here cr x is a Pauli matrix. The cr x ® h operation is 
equivalent to two X operations: X l02) and X (13) , so a CDNOT 
gate is equivalent to two GCX gate. The ,:^ asy of this synthesis 
is 1.60, still less than that in Tab. [Ill]for ququart 1.95. From 
discussion above, it can be seen that our algorithm given here 
is first efficient multivalued synthesis routine which produces 
best known results for all dimensionality d, and for both the 
small n case and the asymptotic case. 

The syntheses of generic n-qubit circuits based on QSD and 
their asymptotic efficiency indexes are listed in Tab. [VI] The 
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TABLE IV: The CINC gate count for the synthesis of qudit quantum circuits obtained by using CSD with balanced partition [ 28h. 





d 




n 

3 

4 

5 

6 

7 

8 

2 

36 

72 

280 

420 

588 

784 

3 

3360 

4464 

69720 

60960 

381780 

142240 

4 

20088 

40824 

670320 

1563660 

4928616 

4750256 

5 

1382952 

685440 

252347440 

38074200 

1.0 x 10 10 

3.1 x 10 8 

6 

8254764 

22254984 

2.5 x 10 9 

3.7 x 10 9 

1.4 x 10“ 

3.9 x 10'° 

7 

127837404 

357389712 

1.4 x 10 11 

1.8 x 10 11 

1.1 x 10 13 

2.5 x 10 12 

8 465572880 

2.9 x 10 9 

8.8 x 10 11 

4.2 x 10 12 

9.8 x 10 13 

8.0 x 10 13 


TABLE V: The CDNOT gate count for the synthesis of ququart quan¬ 
tum circuits & 

n 2 3 4 5 6 n 

Gate count 60 1200 20160 326400 5.2 x 10 6 (5/16)4 2 " — (5/4)4" 


timal for all dimensionality d, so that we can build efficient 
quantum circuits for the qudit d is not a power of two as that d 
is a power of two. It is of great significance to make full use of 
various resources. The algorithm given here is the most effi¬ 
cient qudit synthesis routine so far which produces best known 
results in all respects. 


qubit counterpart of the optimal synthesis for the generic n- 
qudit circuits given here is the qubit circuits based on QSD 
with recursion bottom out at the one-qubit circuit (/ = 1) and 
the optimization for uK n -\(R y ) rotations (the second line of 
Tab. 1 VII) . its is 1.50, less than all asymptotic efficiency 
indexes in Tab. m For the qubit case, there is a most effi¬ 
cient synthesis for generic two-qubit circuits which reaches 
its theoretical lower bound of complexity three CNOT gates. 
The best result for the synthesis of generic u-qubit circuits is 
based on QSD with recursion bottom out at the two-qubit cir¬ 
cuit (l = 2) and two additional optimizations (the fourth line 
of Tab. ED, its ./Z asy is 2.09. Now the asymptotic efficiency 
indexes of the generic n-qudit circuits with d > 5 have been 
greater than this value. Moreover the generic two-ququart cir¬ 
cuit has been more efficient than the most efficient generic 
four-qubit circuit. 

TABLE VI: The CNOT counts of n-qubit quantum circuits based on 
QSD. 

2~~3 4 5 6 n 

1=1 [19] 6 36 168 720 2976 (3/4)x4 2 "-(3/2)x4" 1.33 

/=1, optimal 5 31 147 635 2635 (2/3)x4 2 "-(3/2)x4" + l/3 1.50 

1=2 [19] 3 24 120 528 2208 (9/16)x4 2,, -(3/2)x4" 1.78 

/=2, optimal [19] 3 20 100 444 1868 (23/48)x4 2 "-(3/2)x4"+4/3 2.09 


The multivalued quantum circuits do have advantages over 
the binary quantum circuits. The generic n-qudit circuits with 
d > 5 and generic two-ququart circuits given here are practical 
circuits which are more efficient than the most efficient qubit 
circuits. The leading factor of the lower bound of complexity 
for qudit circuits is small by a factor of d - 1 in comparison 
to that for qubit circuits and the asymptotic efficiency index is 
increased with the increase of dimensionality d, further reveal 
the advantages and benefits of qudit circuits over generic qubit 
circuits. 

There is still plenty of room for improvement in the synthe¬ 
sis of multivalued quantum circuits. One of most important 
work for the improvement is to optimize the two-qudit quan¬ 
tum circuits. Since our algorithm for generic qudit circuits 
given here is recursive, the more efficient generic qudit cir¬ 
cuits can be obtained from more efficient two-qudit circuits. 
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